LABORATORY INVESTIGATION OF DIRECT MEASUREMENT OF ICE 
WATER CONTENT, ICE SURFACE AREA, AND EFFECTIVE RADIUS 
OF ICE CRYSTALS USING A LASER-DIFFRACTION INSTRUMENT 


FINAL REPORT 

December 28, 1995 

NASA AMES RESEARCH CENTER 
(Contract NAS2-14261) 


H. Gerber 1 , P.J. DeMott 2 , D.C. Rogers 2 


1 Gerber Scientific Inc., 1643 Bentana Way, Reston, VA 22090 
2 Colorado State University, Ft. Collins, CO 80523 


LABORATORY INVESTIGATION OF DIRECT MEASUREMENT OF ICE 
WATER CONTENT, ICE SURFACE AREA, AND EFFECTIVE RADIUS 
OF ICE CRYSTALS USING A LASER-DIFFRACTION INSTRUMENT 


FINAL REPORT 

December 28, 1995 

NASA AMES RESEARCH CENTER 
(Contract NAS2-14261) 


H. Gerber 1 , P.J. DeMott 2 , D.C. Rogers 2 


1 Gerber Scientific Inc., 1643 Bentana Way, Reston, VA 22090 
2 Colorado State University, Ft. Collins, CO 80523 




Table of Contents 

Pg- 

— 

Title Page 

1 

— 

Table of Contents 

2 


Summary 

3 

— 

1. Introduction 

4 


2.Theory 

5 


3. Experiment 

8 


3.1 Instrumentation 

8 


3.2 Experimental Procedure 

10 

— 

3.3 Measurements 

11 

* 

4. Analysis 

17 


4.1 Strategy 

17 

u 

4.2 Ice-Crystal Orientation 

18 


4.3 IWC and PSA of Ice Crystals 

19 

ia 

4.4 Effective Radius of Ice Crystals 

21 

i 

4.5 Resolution Limit 

23 


5. Conclusions 

24 


6. References 

25 


7. Acknowledgements 

26 

H 

B 

Appendix A - Experimental Notes 

Appendix B - CSU Dynamic Cloud Chamber Measurements 

27 

~~~ 

for all Experiments 

35 

ma 

Appendix C - Raw Data Produced by the PVM for all 


H 

Experiments 

54 

„ 

Appendix D - Analyzed Data 

73 

5 

Form 298 

142 


2 



SUMMARY 


The aircraft microphysics probe, PVM-100A, designed for measuring liquid 
water content (LWC), droplet surface area (PSA), and droplet effective radius (Re) of 
liquid water clouds, was tested in the Colorado State University dynamic cloud 
chamber to establish its ability to measure related parameters consisting of ice water 
content (IWC), PSA, and Re in ice clouds. The PVM response was compared to other 
means of measuring those ice-cloud parameters in the chamber that included using 
FSSP-100 and 230-X 1-D (Particle Measuring Systems Inc.) optical probes for ice- 
crystal concentrations, a film-loop microscope for ice-crystal habits and dimensions, 
and an in-situ microscope for determining ice-crystal orientation. Intercomparisons 
were made for ice clouds generated in the CSU chamber containing ice crystals 
ranging in size from about 10 um to 150 um diameter, and ice crystals with plate, 
columnar, dendritic, and spherical shapes. 

It was not possible to determine conclusively from the intercomparisons that the 
PVM accurately measures IWC, PSA, and Re of ice crystals, because heat from the 
PVM evaporated in part the crystals in its vicinity in the chamber thus affecting its 
measurements. Similarities in the operating principle of the FSSP and PVM, and a 
comparison between Re measured by both instruments, suggest, however, that the 
PVM can make those measurements. The resolution limit of the PVM for IWC 
measurements was found to be on the order of 0.001 g/m . Algorithms for correcting 
IWC measured by FSSP and PVM were developed. 



1. Introduction 


The purpose of this work is to quantify the operation of the Particulate Volume 
Monitor, Model PVM-100A, in clouds that consist of ice crystals. While this probe 
(termed PVM in the following), designed for aircraft use, was shown to produce 
accurate measurements of liquid water content (LWC), particle surface area (PSA), 
and effective radius (Re) in clouds containing liquid water and droplets with sizes 
ranging over about 3 um to 50 um diameter (Gerber, 1993; Gerber et al., 1994), its 
response to particles of irregular shape, such as ice crystals, is unknown. If such a 
quantification for ice crystals is possible, then the PVM could become a useful tool for 
studying the smaller ice crystals in cirrus clouds and con-trails, because it has 
measurement advantages including a large optical sample volume, and data sampling 
rates of up to 2 KHz. 

The PVM can be classed as a "laser diffraction instrument' 1 , in that it uses a 
collimated laser beam that irradiates suspended particles that scatter light into a 
narrow forward cone. The scattered light consists primarily of light diffracted by the 
particles, as can be shown for the PVM and spherical particles by comparing 
predictions of Fraunhofer diffraction theory and Mie theory (Gerber, 1991). There is 
ample evidence that spherical particles as well as particles of irregular shape diffract 
light proportionally to their projected geometric area (Hodkinson and Greenleaves, 
1983; Hodkinson, 1966; Born and Wolf, 1980; Bohren and Huffman, 1983). A 
motivation for the present work is that the PVM may not respond in that fashion for 
non-spherical particles, because non-linear spatial filters is used in the PVM to weight 
the forward scattered light in such a manner that particle area and volume are 
obtained. This filter may cause the expected proportionality of the diffracted light to the 
projected area of irregular particles to be partially lost. Thus in this work we test the 
following hypothesis: 1) The PVM measures the projected area, the projected area to 
the 3/2 power (volume) of ice crystals, and Re of ice crystals. A second hypothesis is 
also tested: 2) The PVM can estimate values of ice water content (IWC) as small as 
about 0.001 g/m 3 , which is an amount found in diffuse atmospheric ice clouds (e.g., 
Liou, 1992.) 

The method for testing the behavior of the PVM is to expose it to ice crystals of 
known shape and size. This is done in the Dynamic Cloud Chamber (DCC) located at 
Colorado State University in Ft. Collins, CO. This chamber has the capability of 
forming ice clouds similar to those found in the atmosphere, is large enough to contain 
the PVM, and includes a variety of other instrumentation used to characterize the ice 
crystals (DeMott, 1988, 1990; DeMott and Rogers, 1990). Given that the sub-freezing 
temperature and ice-nuclei content can be controlled in the DCC, means that crystals 
with different sizes and habits can be generated (e.g., see Magono and Lee, 1966.) 
The following sections include a theoretical discussion relating diffracted light of ice 
crystals to IWC, PSA, and Re; a description of the experiments run with the PVM in 
the DCC, and data collected; an analysis of the data; and a final section that draws 
conclusions on the validity of the two hypotheses. 
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2. Theory 


* 

We start with the definitions of LWC, PSA, 
water droplets of radius r (see Fig.1): 

and Re for a cloud with spherical 


LWC{glm 3 ) =-^-£ r 3 n{r) 

3 r 

(1) 

— 

PSA{cm 2 lm 3 ) = 4n Y, t 2 r ) 

r 

(2) 

i in a 

. V™ , LWC 

^ p SA 

(3) 


here p w is the droplet density, and n(r) is the droplet size distribution. 




Fig. 1 - Spherical droplet and hexagonal 
ice crystal. 


In developing expressions analogous to Eqs. (1)-(3) for ice crystals, we use as 
an example a solid and symmetrical hexagonal ice crystal shown in Fig. 1, for which 
the length of a facet on the basal plane is a, and the thickness is c. It is obvious that 
the determination of IWC, PSA, and Re for such an ice crystal with a laser-diffraction- 
type optical instrument must include a discussion about the orientation in space of the 
crystal with respect to the axis of the instrument's laser beam. We know from 
atmospheric measurements that ice crystals are sometimes oriented (Platt et al., 
1978). However, it is though that in most situations the crystals are not strongly 
oriented; although, the degree of orientation is usually an unknown. We make the 
assumption here that the ice crystals are randomly oriented in space, so that the 
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crystals will also be randomly oriented with respect to the laser beam in the PVM. 
Evidence of ice crystal orientation for crystals generated in the DCC is given in 
Section 4.2. 

The volume V of the crystal in Fig. 1 is given by (Takano and Liou, 1995) as 

V=Zfea z cf2 (4) 

which can be converted to IWC my multiplication with the density of ice (py — 0.917 
g/cm 3 ). 

They also give 

G= 3a 2 [\/3 + 4{c/2a)}/ 4 (5) 

as the geometric projected area of the randomly oriented hexagonal crystal, where the 
total surface area, PSA, of the crystal is given by 4G. 

In keeping with the form of Eq. (3), a sphere-equivalent Re for hexagonal ice 
crystals, such as in Fig. 1 , is defined by 

Y IWC Y ( \/3 ad2 } n ( % ^ 

Re - 3 p 2.75 (6) 

Y PSA Y {fia+2c) n{a,c) 

a,c a,c 


where the summations must be done over the hexagonal ice-crystal size distribution, 
n(a,c), that depends on both a and c. This definition of Re is consistent with the use of 
Re from the standpoint that the denominator is proportional to the total area 
concentration of randomly oriented crystals affecting the scattering properties of the 
incident radiation, and that the numerator is proportional to the total ice volume 
concentration affecting the absorption of radiation. It would be expected that isometric 
solid hexagon crystals will produce values of Re that are most consistent with Re for 
spheres of the equal volume and area, while hexagonal crystals with significantly 
different values of a and c will differ from Re for spheres, because of differences in 
their scattering cross sections. Small crystals with characteristic dimensions close to or 
smaller than the first Mie peak would show the largest differences. 


For monodisperse hexagonal ice crystals Eq. (6) is simply 


Re- 1.5 & SS- 
s[3a* 2c 


(7) 
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If the PVM operates properly as a diffraction-only instrument it should be able 
to measure G in Eq. (5) and the denominator of Re in Eq. (6) correctly for randomly 
oriented crystals (barring non-linear filter effects.) However, the same is not true for 
the ice volume V in Eq. (4) and the numerator in Eq. (6), because a diffraction 
instrument that is capable of measuring volume a r 3 for a sphere, measures volume 
<* G^ 2 for a randomly oriented particle of irregular shape. The latter diffraction volume, 
V d , is defined by 


1/,= 47T 1/2 6 3/2 /3 


(8) 


V d is larger than the true value of the ice volume V of the crystal, so that the 
PVM will overestimate IWC as well as the sphere-equivalent Re of the hexagonal ice 
crystals. It is thus necessary to apply a correction factor to the IWC and Re measured 
by the PVM. This factor, C, is given for solid hexagonal crystals by 

q_ V_ _ 3\/n_{c/d) (g) 

Kr [3^3 + 2 (cidjf 2 

The value of C given by Eq. (9) is shown in Fig. 2 as a function of the dimensions c 
and a of hexagonal as well as other shapes of randomly oriented ice crystals. Figure 2 
suggests that for a range of c/a typical of atmospheric hexagonal ice crystals the PVM 
should overestimate IWC and Re by a factor of between about 1 .25 to 2. 



Fig. 2 - Volume correction factor for a diffraction instrument for 
randomly oriented ice crystals. Details given in text. 
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The foregoing illustrates a set of equations developed for IWC, PSA, and Re for 
solid, symmetrical, and randomly oriented hexagonal ice crystals measured by a 
diffraction instrument. Atmospheric ice crystals are often more complex than the 
simple crystals discussed so far. Such complex crystals can include dendrites, capped 
columns, and bullet rosettes, which can also have various degrees of hollowness. 
Takano and Liou (1995) give expressions for V and S for models of these more 
complex randomly oriented crystals, so that equations similar to the preceding can be 
derived for IWC, PSA, and Re. Calculated values of C for these other crystals are 
shown in Fig. 2. The parameters bt/a for dendrites, a3/a2 for capped columns, and d/c 
for hollow columns shown in Fig. 2 are defined in Takano and Liou (1995); the 
"shadowing factor", s, is applied to all three crystal types; and the crystal dimension 
ratio c/a for capped columns is defined here as (l 1 +l 2 +l 3 )/a. The numbered data points 
are described in Section 4.4. 

This analysis shows that, if the PVM is not adversely affected by the non- 
linearity of its spatial filters, it should overestimate IWC as well as Re, while giving the 
proper value of the projected geometric area of ice crystals. The work by Takano and 
Liou (1995) provides the mechanism by which these overestimates can be corrected, 
if the habit and orientation of the ice crystals are known. 


3. Experiment 
3. 1 1nstrumentation 

The instrumentation used for testing the PVM-100A is located at the Colorado 
State University Cloud and Aerosol Simulation Laboratory. AD tests are conducted 
using the Dynamic Cloud Chamber (DCC) located at that laboratory, with most other 
instruments placed on or inside the DCC. The DCC and equipment used in the tests 
are summarized in Table 1. 


Table 1 - Instrumentation used for testing at the CSU Cloud and Aerosol Simulation Laboratory the 
response of the PVM to ice crystals. 


Instrument 


Description 


1. Dynamic Cloud Chamber 


1.2 m 3 pressure vessel 1 
cooling capabilities, wi 
for probes and aerosol 
Range of Operation: 
Temperature 
Pressure 
RH 

Vertical Velocity 


ith air-expanding and 
i vents and windows 
nd gas injection. 

+40C to -55C 
90 kpa to 50 kpa 
0.1% to >100% 

0.2 m/s to 20 m/s 
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2. Thermocouples (Type-E) 

3. Strain Gauge 

4. G.E. and E.G.&G. hygrometer 

5. TSI Nucleus Counter (3020) 

6. Film-Loop Microscope 

7. In-Situ Microscope 

8. FSSP-100 

9. 230-X Probe 

10. PVM-100A 


11. Laser 


inside wall and air temperature of DCC 

pressure inside DCC 

optical dew-point hygrometers for RH 

measurement in DCC 

ice and condensation nucleus 

concentration 

continuous horizontal 16 mm film loop near 
bottom of DCC on which ice crystals 
sediment and are observed with a 
microscope/CCD video and recording system 
long-focal length microscope (Infinity, K2) for in-situ 
observation of crystals in DCC, observed with 
CCD video and recording system 
Forward-Scattering Spectrometer probe (PMS), 
particle spectra from 2 um - 47 urn diameter 
1-D Spectrometer (PMS), particle spectra from 
20 um - 300 um diameter 
Particulate Volume Monitor (GSI), LWC, PSA, 

Re; size sensitivity 3 um - 50 um (100% points), 

2 um - 70 um (50% points) 

diode laser (680 nm) used to observe backscatter and 
scintillation from crystals 


The locations of the optical instruments in the cylindrical DCC are as follows: 
the film-loop microscope is near the base plate of the DCC; the FSSP-100 is located 
about 1/5 up the outside wall of the DCC, sample air is drawn from the chamber at a 
low rate; the 230-X is located below the bottom plate of the chamber, sample air is 
drawn from the chamber; the in-situ microspcope is mounted in and about 1/5 up the 
chamber wall; and the PVM is located for most tests in the chamber so that the focal 
point of the in-situ microscope is positioned just above the cylindrical aperature of the 
PVM, see Fig. 3. The PVM is further positioned so that the opening of its cylindrical 
aperture and sensitive volume is horizontal to permit unobstructed sedimentation of 
the crystals. For some of the tests the PVM is moved more towards the center of the 
DCC. The exposed metal shell of PVM is thermally insulated with light-weight closed- 
cell styro-foam in order that heat exchange between it and the chamber environment 
is kept at a minimum. 
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Fig. 3 - Location of PVM inside the Dynamic Cloud Chamber 
(DCC) with respect to the focal point F of the in-situ microscope. 


3.2 Experimental Procedure 

The goal for generating ice crystals in the DCC is to create ideally 
monodisperse ice crystals varying in size and habit for different experiments, in order 
that the evaluation procedures for the optical instrument would be less time 
consuming. It is desired to generate crystals ranging in size from about 10-um to 200- 
um diameter, and, at the minimum, to generate crystal habits consisting of plates, 
dendrites, columns, and spheres. The experimental procedure used in attempting to 
meet this goal is summarized as follows: 

The chamber pressure and temperature is reduced by computer control at a 
constant rate, simulating a moist adiabatic ascent of about 5 m/s. Given that the air in 
the DCC originally contains room moisture, a cloud forms when the dew point is 
reached. Upon reaching a target sub-freezing temperature, chosen from 
temperature/ice-crystal habit relationships (e.g., Magono and Lee, 1966), the 
temperature of the supercooled cloud is kept constant. Ice nuclei of a given 
concentration are then injected into the chamber to nucleate water vapor and 
supercooled water to form a totally glaciated cloud consisting only of ice crystals. The 
size of the resulting ice crystals is roughly proportional to the number of injected ice 
nuclei, as well as the dew point of the air prior to condensation. An attempt is made to 
operate all instruments in the chamber during each entire experiment; however, a 
detailed analysis of the data is only done here for intervals of time following total 
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glaciation of the cloud. 


3.3 Measurements 

Eighteen experiments in the DCC were run over a period of two weeks. The 
experiment # is identified as PVM#, and each interval #, lasting 60 sec, for which a 
detailed analysis is done on the glaciated cloud is given by I#. The experiments are 
summarized in Table 2. 


Table 2 - Summary of experiments run in the CSU Dynamic Cloud Chamber 


Exp. No. 
PVM# 

Interval No. 
I# 

Detailed Ice- 
Crystal Analysis 

Comments 

1 


no 

no seeding 

2 


no 

no film loop 

3 


no 

marginal LWC and xls 

4 


no 

marginal LWC and xls 

5 

1,2,3 

yes 

hex plates with some dendritic structure, 
up to about 60-um diameter. -10C 

6 

4 

yes 

hex plates to about 150-um, -10C 

7 


no 

no film loop 

8 


no 

no film loop, stobe and in-situ microscope 
properly operating 

9 

5, 6 

yes 

hex plates to about 50 um, -18 C; add 
aspirator to base of PVM 

10 

7 

yes 

columns to about 50 um , -6 C 

11 

8,9 

yes 

mostly hex plates < 20-um, -23.5 C 

12 

10, 11 

yes 

hex plates with dendritic features, to 100 um, 
-15.5 C 

13 


no 

overseeded, columns and plates, -10C 

14 


no 

no film loop 

15 

12 

yes 

hex plates to 80 um, -12 C; move PVM 
towards center of DCC 

16 

13 

yes 

columns to 40 um, -7 C 

17 

14, 15 

yes 

spherical poly-crystals about 10 um, -30 C 

18 

16, 17, 18 

yes 

mostly columns, - 33 C 


Detailed notes describing each experiment are given in Appendix A. The raw 
time-dependent data consisting of temperature, pressure, RH, and FSSP and 1-D 
probe outputs measured in the DCC are illustrated for one experiment in Fig. 4; 
Appendix B shows this data for all experiments. Figure 5 shows the raw data output 
by the PVM for the same experiment; Appendix C shows the PVM ouput for all the 
experiments. 


li 
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0 500 1000 

USM (mc) 


Fig. 4 - Pressure, temperature, RH, and FSSP and l-D probe 
data measured as a function of time during Experiment 
PVM5. 
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Fig. 5 - PVM measurements of IWC, PSA, and Re during Exp. PVM 5. 
The scale of the vertical axis gives IWC (g/m 3 ), Re/10 (urn), 
and is relative for PSA. Intervals 11-13 for which detailed ice crystal 
analyses are done are shown. 
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Figure 6 shows examples of images of ice crystals collected with the film-loop 
microscope and the in-situ microscope. Continuous video recordings by both systems 
were made during all tests. The distribution of a and c dimensions of crystals were 
determined by manually measuring large numbers of images on the film-loop video. 
The in-situ microscope images were used to determine the orientation of the crystals 
as they sediment just above the open cavity of the PVM. Crystals of single habits 
were generated as expected; however, they were not monodisperse. 




Figure 6 - Examples of video images collected with the film-loop 
microscope (A, Exp. PVM9; B, PVM 16), and with the in-situ 
microscope (C, PVM7; D, PVM8; E, PVM8). The grid spacing 
in A, B is 200 urn; the scale in C,D,E is estimated from the 
50-um crystal diameter in C. 
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The raw data for the optical probes, FSSP, 1-D, and PVM, are reduced and 
assigned new variable names as defined in Table 3. Figures 7-13 illustrate the 
nature of the reduced data for experiment PVM5. Appendix D gives similar sets of 
graphical data for all the experiments for which a detailed analysis was done on the 


Table 3 - Definitions of variables used in the description of the reduced data shown in 
Figs. 6-12, and in Appendix D. 

P VMWC: PVM 1 00-A 1 0 Hz liquid or ice water content. 

PVMWC5s: Running 5s average of PVMWC. 

TOTWC: Total calculated liquid and/or ice water content. 

FVIDIWC. FSSP ice water content computed by combining the video film loop size 
distribution, mass versus size relations from Redder and Fukuta (1989), 
and 5s average FSSP concentration. 

1DVIDIWC: 23 OX (ID probe) ice water content computed by combining the 
video film loop size distribution, mass versus size relations and 5 s 
average ID concentration 

GELWC: liquid water content estimated from the heated-inlet dewpoint 

hygrometer. 

FRFIWC: FSSP ice water content computed using the empirical formulas of Redder 

and Fukuta and 5s average FSSP size distribution to represent major axis 
dimension. 

FCORIWC: FSSP ice water content computed as for FRFIWC, but outside the 

regime of the Redder and Fukuta formulas. Ice volume was based on 
the observed ice crystal shape and an assumed ice density pi = 0.7 g m 3 . 

1DCORIWC: ID ice water content computed outside the regime of the Redder and 

Fukuta formulas, based on explicit shape of ice crystals and an assumed ice 
density pi = 0.7 g m' 3 . 

PVMPSA: PVM 1 00-A particle surface area (PSA) channel. 

PVMPSA5s: Running 5s average of PVMPSA. 

TOTPSA: Total calculated PSA. 

FVIDPSA: FSSP PSA calculated based on the video film loop size distribution 

1 DVIDPS A: 1 D PS A calculated based on the video film loop size distribution. 

FRFPSA: FSSP PSA calculated using the FSSP size distribution to represent major 

axis dimension and Redder and Fukuta axial relationships. 

1DRFPSA: ID PSA calculated using the ID size distribution to represent major axis 

dimension and Redder and Fukuta axial relationships. 

PVMREFF: PVM1 00-A effective radius. 

FSPHREFF: FSSP effective radius calculated assuming the particles are measured as 

spheres. 

1DSPHREFF: ID effective radius calculated assuming the particles are measured as 
spheres. 

TOTREFF; Total calculated effective radius using TOTWC and TOTPSA. 

I##: 60 s interval number (##) of particular interest during which only ice 

crystals were present and video size distributions were measured. 
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glaciated cloud in the DCC. All sets, except for experiments 11, 17, 18 contain the 
relative distributions of a and c for the crystals measured manually from the film loop; 
the exceptions are due to the small size of the crystals which sediment inefficiently 
onto the film. In these three cases, FSSP data was used for establishing ice crystal 
sizes. The reduced-data variables in Figs. 7-13 are defined and discussed further in 
the next section. 


PVM5 

Water Content Comparison 



Fig. 7 - Comparison of 
LWC and IWC cal- 
culated and measured 
by FSSP, 1-D, and 
PVM for exp. PVM5. 


PVM5 

PARTICLE SURFACE AREA 



Fig. 8 - Comparison of 
measured and cal- 
culated PSA. 
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FRACTION Cone, (cm- 1 urn' 1 ) Reff(um) 


PVM5 

EFFECTIVE RADIUS 



PVM5 

FSSP vs. 1-D Size Distributions 



Fig. 9 - Comparison of Re 
measured with FSSP, 1-D 
and PVM. 


Fig. 10 - FSSP and 1-D 
ice crystal size distrib- 
utions. 


PVM5 

Comparitive Fractional Size Distributions, Interval 1 



Fig. 1 1 - Relative size 
distributions of a and c 
crystal axes determined 
from film-loop video, 
FSSP and 1-D. 
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PVM5 

Video Axis Relationship 



Fig. 12 - Film-loop 
video measurements of 
crystal dimensions a 
and c, and the 
predictions of the 
Redder-Fukuta (1989) 
realtionship for two 
temperatures. 
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Fig. 13 - The ratio of 
total LWC and IWC to 
PVM LWC and IWC. 


4. Analysis 
4. 1 Strategy 

The requirement for comparing PVM outputs to measured values in the DCC of 
IWC, PSA, and Re for ice crystals can be addressed in various ways, given that four 
other optical systems, FSSP, 230-X, in-situ microscope, and film-loop microscope, are 
used in the DCC to collect size and concentration information on the ice crystals. 

None of the systems is judged to be able to produce by itself accurate values of IWC, 
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etc., for comparison with the PVM: The FSSP and 230-X probes respond to the non- 
spherical ice crystals in an unknown fashion, the film-loop microscope gives only 
relative values of ice-crystal concentration, and the very small sample volume of the 
in-situ microscope results in insufficient data for concentration measurements. 

It was decided that the most suitable role of the instruments for determining 
IWC, PSA, and Re of the crystals was to use the film-loop microscope to establish, by 
manual measurements of the ice crystal dimensions on the film-loop video, the relative 
size distribution of the ice crystals, and to use the total ice crystal concentrations 
measured by FSSP and 230-X probes to make these relative distributions absolute. If 
the crystals were smaller than about 42-um diameter, then only the FSSP 
concentrations was applied, and when the crystals size distributions extended above 
that size, the FSSP and the 230-X concentrations were both applied. This procedure 
is illustrated in Fig. 1 1 where the relative size distribution of the crystal dimensions a 
and c, as measured on the film-loop video, is shown for Experiment PVM5; and Fig. 

10 illustrates the FSSP and 230-X concentrations that are applied to the relative 
spectra in Fig. 11. 

The procedure to obtain IWC from the ice crystal sizedistributions is to apply 
empirical algorithms (Redder and Fukuta, 1989) that relate ice crystal dimensions to 
ice crystal mass. This is done for all experiments except for PVM11, PVM17, and 
PCM 18 for which the crystals were too small to efficiently sediment onto the film loop 
(see Table 3.) For those three experiments the FSSP size distribution data is used 
instead; this is felt to give reasonable results, because the crystals in those three 
experiments are either nearly isometric or spherical. IWC values determined in the 
preceding ways are labeled TOTWC in the reduced data (see Fig. 7, and Appendix D.) 

PSA is also determined by relying on the relative size distribution of the ice 
crystal dimensions measured on the film loop multiplied by total FSSP and 230-X 
concentrations. PSA is calculated by applying the relationship for PSA = 4 x G given 
by Takano and Liou (1995) for randomly oriented ice crystals, see Eq. (5). PSA values 
determined in this manner are labeled PSATOT in the reduced data (e.g., see Fig. 8.) 

Re is determined by using Eq. (6) with TOTWC and TOTPSA. Re values 
determined in this fashion are labeled TOTREFF in the reduced data (e.g., see Fig. 9.) 

4.2 Ice-Crystal Orientation 

The treatment in the preceding section requires the assumption that ice crystals 
produced in the DCC are oriented randomly in space. While an absolute statement 
can not be made about the orientation of all the crystals, the following evidence is 
presented that suggests random orientation: 1) The in-situ microscope showed ice 
crystals sedimenting with their longer axis at various angles with respect to the 
horizontal (e.g., see Fig. 6.), 2) ice crystals sedimented onto the film-loop with their 
axes at various angles with respect to the horizontal surface of the film (this includes 
plates some of which land on edge that permits a measure of the c dimension), and 3) 
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the back-scattered light from the diode laser projected into the ice cloud showed 
strong scintillation suggesting the presence of tumbling ice crystals. 

4.3 IWC and PSA of fee Crystals 

Analyzed data consisting of measured values of IWC and PSA are 
summarized in Table 4 (see Table 3 for definition of variable names.) Table 4 shows 
that not all experiments produced analyzed data, as also indicated in Table 2. The 
criteria for choosing to analyze an experiment's data depended largely on the proper 
operation of the film-loop microscope, which did not function occasionally. 

Overseeding or underseeding with ice nuclei also caused non-analysis of the 
experiment. 

A typical result in comparing IWC (defined as PVMWC) measured by the PVM 
with TOTWC (see Section 4.1) is shown in Figs. 7 and 13 for Experiment PVM5. 
PVMWC is substantially less than TOTWC. The maximum TOTWC agrees reasonably 
well in Fig. 7 with the theoretical IWC calculated by differencing the appropriate vapor 
mixing ratios in the chamber. This suggests that the PVM was underestimating IWC 
for some reason; whereas, it should have been overestimating IWC according to the 
arguments presented in the theory section. The underestimate is also evident in the 
liquid cloud as Fig. 13 shows. Figure 14 summarizes the comparison between 
PVMWC and TOTWC. Except for some columnar crystals of low IWC, the trend in the 
underestimate of IWC by the PVM is consistent for all experiments. 

We concluded that the PVM underestimated both IWC and LWC during the 
experiments, because we were unable to eliminate entirely heat conduction from the 
PVM during the experiments. This caused water droplets as well as ice crystals to 
evaporate in the vicinity of the PVM as the chamber was cooled. The diode laser 
confirmed this behavior, because the backscattered light from ice crystal clouds 
showed a minimum area located just above the circular opening of the PVM during the 
experiments. This evaporative effect was not complete, as shown by the existence of 
some backscatter just above the PVM, as well as by the outputs from the PVM; 
however, the exact nature of the evaporation remained unknown. It appeared that the 
cylindrical annulus of the PVM acted as a warm chimney, despite careful insulation of 
the metal PVM probe with styrofoam. This effect became apparent early in the course 
of the two-week experiment. In an attempt to counteract the heating in the annulus a 
low-speed suction fan was constructed and attached to the bottom opening of the 
annulus to draw the crystals more rapidly past the PVM heat sources. This was done 
starting experiment PVM9. The results were largely unchanged. A final attempt to deal 
with this problen was done starting experiment PVM15 when the PVM was moved 
closer to the middle of the DCC, and it was raised above the chamber floor to a height 
of 9". The latter modification gave the best comparisons between PVMWC and 
TOTWC, especially for experiments PVM 17 and PVM 18. This improvement may, 
however, have been due to the slowed evaporation, because of the low ice vapor 
pressures existing at the cold temperatures used in those experiments. 
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Fig. 14 - Comparison of IWC measured by the PVM (PVMWC) 
and determined from film-loop crystal measurements 
(TOTWC, see Section 4.1.) for the experiments listed 
in Table 4. 


The same evaporative effect caused by the PVM is evident in the comparisons 
between PVM PSA and TOTPSA, as illustrated in Fig. 8. The evaporation could be 
related to inhomogenous or homogenous mixing, or a combination of the two. In the 
event that the first possibility holds, then the value of Re measured by the PVM may 
have remained unchanged, as will be discussed in the next section. 

4.4 Effective Radius of ice Crystals 

Table 4 lists the four measurements of effective radius, Re, of the ice crystals 
that were made in the cloud chamber during the experiments: FSHREFF (FSSP-100), 
IDSPHREFF (230-X), PVMREFF (PVM-100A), and TOTREFF (film-loop microscope 
and total ice crystal concentrations from PMS probes.) Re for the FSSP and the 230-X 
is calculated by assuming spherical particle shape. 

Figure 15 compares Re measured with all the above probes except the 230-X, 
which gave Re values that were obviously too large, because most of the ice crystals 
had sizes smaller than the 42-um lower detection limit assigned to this probe. The top 
part of Fig. 15 shows the data as listed in Table 4, with the PVM giving apparent 
excellent agreement with the TOTREFF values, and the FSSP giving values that are 
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TOTREFF (um) 



TOTREFF (um) 


Fig. 15 - Comparison of effective radius (REFF) measurements of 

ice crystals made with PVM and FSSP vs film-loop micros- 
cope measurements. The corrected plot applies Eq. (9). 

mostly too large. The bottom "corrected" part of Fig. 15 gives the Re values for the 
PVM and FSSP listed in Table 4 adjusted by the correction factor C given by Eq. (9) 
for which average measured values of c and a are used for all experiments except 
PVM 12. The crystals were assumed to be solid hexagons, and the shadowing factor 
was not applied. These values of C are shown in Fig. 2. Applying C results in the PVM 
underestimating TOTREFF in the bottom part of Fig. 15, while the FSSP shows good 
agreement on the average with TOTREFF, as indicated by the average ratio of 
FSHPREFF/TOTREFF = 0.98. The spread in the data point may be less if the 
contributions to the actual values of Re from the larger crystals measured by the 230- 
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X, present in some of the experiments, are included. 

We conclude that the good agreement between PVMREFF and TOTREFF in 
the top part of Fig. 15 is fortuitous, and that the corrected portion of Fig. 15 shows the 
more likely relationship between these Re measurements. This result further suggests 
that the FSSP appears to behave primarily as a "laser diffraction" instrument for which 
corrections for randomly oriented ice crystals similar to those given by Eq. (9) can be 
applied. The evaporation effect caused by the PVM must have been at least in part 
related to a homogenous mixing process, because the values of Re were not 
conserved as shown by the lower portion of Fig. 15. 

4.5 Resolution Limit 

The resolution limit of the PVM for measuring IWC was determined by 
observing the performance of the instrument during a period in an experiment when 
the ice crystal concentration was very small. Near the end of Experiment PVM 18 such 
a period occured and is show in Fig. 16. 



Fig. 16 - IWC measured by the PVM during Experiment PVM18. 


Figure 16 shows that the PVM is capable of detecting IWC values on the order 
of 0.001 g/m 3 . The spikes indicate the passage through the instrument's laser beam of 
individual ice crystals. The nature of the operating principle of the PVM is such that 
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individual spikes do not indicate accurate IWC values; however, a sufficiently long 
average of the spikes over time results in meaningful IWC values. 

5. Conclusions 

1) We were unable to demonstrate conclusively that the non-linear spatial filters 
used in the PVM to weight scattered light do not prevent the PVM from measuring the 
projected area, the projected area to the 2/3 power, and the effective radius Re of ice 
crystals generated in the CSU cloud chamber, because the residual heat conducted 
from the probe evaporated a portion of the crystals near it, and affected the 
measurements. However, the FSSP measurements made on the same crystals, where 
heating was not a problem, suggest indirectly that the PVM should be able to make 
those measurements. 

2) After correcting the FSSP measurements to convert "diffraction volume" to 
ice crystal volume, this probe gave values of Re in best agreement with the Re 
established by manual sizing of the ice crystal various shapes and sizes sedimented 
onto the film-loop microscope located in the chamber. This suggests that FSSP 
measurements are useful for Re measurements when crystal sizes fall into its range of 
operation, and that scattered-light dependencies as a function of scattering angle other 
than those proportional to crystal size 2 and size 3 , which the FSSP and also the PVM 
have, do not affect strongly the estimate of Re. 

3) The measured resolution limit of the PVM for IWC of ice crystals is on the 
order of 0.001 g/m 3 . This limit depends on the noise environment in which the probe is 
used. 


4) A "sphere-equivalent" Re of the ice crystals can be defined for laser- 
diffraction instruments such as the PVM and FSSP when the "diffraction volume" 
measured by the instruments is corrected to give ice crystal volume, and the 
instruments measure the projected surface area of the ice crystals assumed to be 
randomly oriented in space. Correction factors for volume are derived and given for 
typical ice-crystal habits. 

5) The PVM and FSSP are not stand-alone instruments for measuring IWC, 
PSA, and Re of ice crystals, unless approximate values of those parameters are 
desired. These approximate values should be determined using mean correction 
factors described here, and require that the majority of the crystals have sizes that fall 
within the instruments' range of operation that extends to about 50-um diameter. 
Improving the accuracy of those parameters is directly related to the application of 
additional information on the habit, large-size fraction, and orientation of the ambient 
ice crystals. Such information can be measured, for example, with replicators and the 
holographic instrument from Spec. Inc. It is expected that the PVM will produce useful 
and unique data on smaller ice crystals when used in conjunction with those other 
methods, because of its relatively large sampling volume and fast response time. 

When large and small ice crystals are found at the same time in the atmosphere, PVM 
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measurements should be combined with measurements from probes such as the 1-D 

and 2-D instruments from PMS. 
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6/8/95: 

PVM1 


Ascent rate = 5 m s' 1 
Tdew = -2°C 

CCN bubbled 1 minute NaCl 

Mix fan on high speed until later in test, then off 

No video film loop for ice crystal habit 

PMS FSSP and lD-c not on until later in tests because reconnection of interface card was 
not made properly 


6/9/95: 

PVM2 


Ascent rate = 5 m s' 1 
Tdew = 5°C 

CCN bubbled 3minutes NaCl 

Mix fan on low speed 

No video film loop for ice crystal habit 

Seeded with Agl nuclei at approx. 1549(-15°C); cloud rapidly glaciated 
ID had no flow until 1551 (replacement pump not reconnected properly) 

PVM had good response to ice formation 

Orientation video operating. Focus set for 6.5x4.75, 200|.im squares. Video review 
showed that off-forward-axis lighting was causing refraction through ice crystals, so that 
lighting is variable and details are not discernable. We should switch to forward lighting 
again For the crystal sizes nucleated, could also stand a larger FOV. 


6/12/95: 

PVM3 


Focus on orientation video reset for 6.5x8 200>.im squares before camera reinstalled 
Ascent rate = 5 m s' 1 
Tdew = -2°C 

CCN bubbled 3minutes NaCl 
Mix fan on high speed entire test 

Planned to seed with Agl nuclei at approx. -10°C and halt expansion 
Seeding was at 154900, ascent stopped at 154920. 

Either cloud did not glaciate as much as expected, or ice did not grow sufficiently in size 
Video film loop operating this test, but record started late 
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Following the continuous expansion, about 5 short rapid expansions were done to 
nucleate and grow more ice crystals These temporarily cooled the air (not the walls): 

1553 25mb (10+15 in succession) 

155415 10 mb 

155445 10 mb 

155615 10 mb 

155648 10 mb 


6/13/95: 

PVM4 


Focus on orientation video reset for 6.5x8 200pm squares before camera reinstalled. 
Ascent rate = 5 m s' 1 
Tdew = 5.5°C 

CCN bubbled 3minutes NaCl 
Mix fan off entire test 
Both videos operating 

Planned to seed with Agl nuclei (2x dilution of AgI-NH 4 I combustion aerosol) at approx. 
-10°C and halt expansion for plates only 
Seeded at 153620 (-10°C air temperature) 

Stop ascent at 153710 (-12°C wall temperature) 

Dewpointers unbalanced flows (last 3 tests) fixed for this test and heat added to sample 
inlet. EG&G unit looks particularly good for LWC retrieval when cloud was present. 

Test appeared competitively seeded, that is, not entirely glaciated. Excellent plate records 
on video film loop No orientation video because strong lighting washed out images. 


6/14/95: 

PVM5 


Focus on orientation video reset for 6.5x8 200pm squares before camera reinstalled. 
Ascent rate = 5 m s' 1 
Tdew = 5°C 

CCN bubbled 3minutes NaCl 

Mix fan was left on low until just after seeding 

Both videos operating, strobe system readjusted for this test 

Planned to seed with Agl nuclei (2x dilution of Agl-NFLJ combustion aerosol) at approx. 
-10°C and halt expansion for plates only, heavier seeding 
Cloud formation noted at 1 12000 
Seeded at 1 12650 (-10°C air temperature) 
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Orientation video appears to indicate random crystal orientation. Excellent video film loop 
of plate crystals Higher concentrations collected by film loop, in agreement with stronger 
ID signal. 


PVM6 

Ascent rate = 5 m s' 1 
Tdew = 6°C 

CCN bubbled 3minutes NaCl 

No mixing fan 

Both videos operating 

Planned to seed with Agl nuclei (2x dilution of Agl-NHJ combustion aerosol) at approx. 

-10°C and halt expansion for plates only, heavier seeding 

Cloud formation noted at 142720 

Seeded at 1436 (approx. -14°C air temperature) 

Stopped expansion at 143630 

Noted that slight leak through CCN bubbler system allowed some room air IN into 
chamber before seeding was done. 

PVM7 

Ascent rate = 5 m s' 1 
Tdew = 4°C 

CCN bubbled 3 minutes NaCl 

No mixing fan 

Both videos operating? 

Seeding planned for -14°C. 

Cloud formed at -1608, -2°C 
Seeded at 161530, -16°C 
Stopped ascent at 1616 

Video film loop stuck in position at some point. 


6 / 15 / 95 : 

PVM8 

Planned cloud formation test only after strobe system changed to a pinhole light source. 
Ascent rate = 5 m s’ 1 
Tdew = 8 5°C 

CCN bubbled 3minutes NaCl 
No mixing fan after 1030 
Both videos operating 

Dewpoint inlet was unheated, but did not seem to affect EG&G, so enough heating from 
room around sample line. 

Seeded at 104350 on data system computer 
Stopped ascent at 104545, — 18°C 
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8 I 111 


Excellent orientation video shoeing diffraction patterns. 
No video film loop because it stuck in position 


Sri 


6/16/95: 

PVM9 

Ascent rate = 5 m s' 1 
Tdew = 9°C 

CCN bubbled 3 minutes NaCl 

No mixing fan, except on new PVM sampling system. This system consisted of a laminar 
flow element placed below the PVM with a fan for aspiration of the sample volume. 
Should provide 1-2 cm s' 1 aspiration. Voltage on the brushless fan was set to 6V this test 

Both videos operating 

Seeded at 101731, -18°C 

Ascent stopped -101810, T=-19.5°C 

Problems this test: 

1. ID tubes not reconnected initially, so no measurements made early in cloud formed, and 
the pressure control oscillated to make up for the incoming air. A spike in the ID data was 
caused when the 1 D was finally turned on. 

2. The 6V setting on the PVM aspiration fan was clearly not enough. No voltage gave 
LWC=0. Better LWC values were obtained with the voltage set to 12V. I suggest miving 
PVM higher in chamber volume and not using aspiration. 

3. The video film loop was started late and it was not moving at early times during cloud, 
so artifacts were noted in the first few minutes of recorded data 

4. The orientation video was not plugged in until well after cloud formed. 

PVM 10 

Ascent rate = 3 m s‘ l 
Tdew = 9.5°C 

CCN bubbled 3minutes NaCl 
No mixing fan 
PVM sampling fan to 12V 
Both videos operating 

Slight offset on initial temperature profile due to ascent profile memory problems in 
starting ascent 

Seeded at 1 13710, T=-6°C, did not glaciate cloud 
Stopped ascent at -10°C 

Lots of columns on film loop One in focus on orientation film at 1 14045. 
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PVM11 

Attempt to get column/sheath type crystals below -24°C 
Ascent rate = 3 m s' 1 
Tdew - -4°C 

CCN bubbled 3minutes NaCl 
No mixing fan 
PVM sampling fan to 12V 
Both videos operating 

Dewpoint temp rose during early part of expansion, so dry air was not well mixed 
throughout pressure vessel initially. Cloud could be superadiabatic initially. 

Cloud formed at 145747, T = -6.5°C air temp 

Low concentration cloud formed (poor CCN inject?). LWC<0.05 g m' 3 
FSSP switched range 0 to 1 at -150330 

Dewpoint relative humidity appears offset low this test, never registered cloud. 

A few ice crystals present throughout cloud lifetime. 

Video film loop stuck occasionally, more so later in experiment. 

Cloud disappeared by -22°C. 

Seeded anyway at 151025, -23.5°C 
Stopped ascent at T=-27. 6°C 

~17mb pop expansion at 151240 formed all ice cloud centered at 10-1 2ptm 
13 mb more at 1515. Ice grew to 18-20).im 

ID did not.respond well, may have iced over inlet or very small (<30|.im) crystals present. 

Strobe system had problems until aobut 1512 when Brian noticed loose connections 
Some good images both before and after this time. 

Video film loop shows mostly columns, few plates. 


6/19/95: 

PVM12 

Ascent rate = 5 m s’ 1 
Tdew = 6°C 

CCN bubbled 3minutes NaCl 
No mixing fan 

PVM sampling fan to 1 2 V. Not on until 104115 
Both videos operating 

FSSP on late (PMS switch was off) so missed first few 1Hz cloud records 
Video film loop also started recording about 30s late for ice formation. 
Seeded at 104750, T=-15.5°C 
stopped ascent 104850 



PVM13 

Ascent rate = 5 m s' 1 
Tdew = 6°C 

CCN bubbled 3minutes NaCl 

No mixing fan 

PVM sampling fan to 12V. 

Both videos operating 

Seeded at 115435, T=-10°C air, continued cooling to -15°C 
Ice was present early at -8°C (columns) 

Seeding was overseeding. Plates of all sizes and thicknesses apparent. Some were so thick 
that I believe they look like columns. 

Note noise in data system last tow tests must be omitted in plotting data. 

PVM 14 

Ascent rate = 5 m s' 1 
Tdew = 9.5°C 

CCN bubbled 3minutes NaCl 

No mixing fan 

PVM sampling fan to 12V. 

Both videos operating 
Attempt to get larger crystals 

Seeded at 143910 
Stopped ascent at 143940 

Noted crystals formed to >150(.tm on ID updates, with a mode at 70-80pm 
Video film loop stuck early in experiment, so no habit data! 


6/20/95: 

PVM 15 

Repeat last test with 1 additional dilution of aerosol 
Ascent rate = 5 m s' 1 
Tdew = 8.5°C 

CCN bubbled 3minutes NaCl 
Mixing fan on high 

PVM fan removed and PVM raised an additional 5 inches to sit 9 inches off chamber 
floor, very near to the center of the chamber. 

Video film loop only this test while Brian works on strobing faster 

Seeded at 143235 

Film loop stuck some, so some crystal melting noted at times. Otherwise excellent images 
of flat and on-end plates for IVVC calculation Good test for larger crystals. 
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Appears that residua! heating from PVM may be decreasing its LWC value compared to 
the hygrometer and FSSP measurements. 

PVM 16 

Repeat last test, but seed at -7°C for columns 
Ascent rate = 5 m s' 1 
Tdew = 8.5°C 

CCN bubbled 3minutes NaCl 
Mixing fan on high 
video film loop only 

Seeded at 152145 

Excellent test execution. Only columns noted with a final air temperature of -8.5°C 
Good film loop data for both size and concentration. 


6/21/95: 

PVM 17 

Ascent rate = 5 m s' 1 

Tdew = -10°C at start, warmed to -7.5°C during the experiment before falling again 
CCN bubbled 3minutes NaCl 
Mixing fan on high 
video film loop only 

Attempt to seed for columns below -24°C. Lots of ice formed before seeding. Was 
contamination in the air and from the CCN bubbler. 

Seeded at 100320 

Fan was turned off near the end because the small crystals did not appear to be efficiently 
captured on the film loop with the fan on high 

Excellent FSSP/PVM comparison for polycrystalline, nearly spherical, ice crystals. 
Dewpoint hygrometer data not of much use at this colder temperature. Time lag and lack 
of secondary cooling are problems. 

PVM 18 

Ascent rate = 5 m s' 1 

Tdew = -15°C at start, warmed a little during the experiment before falling again 
CCN bubbled 3minutes NaCl 
Mixing fan on low speed until seeding 
video film loop only 

Solid cloud formed at -15°C 
Seeded at 1 1 5430 (— 33°C) 

Mostly columns and some morph spherical particles noted on video. 
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APPENDIX B 


CSU Dynamic Cloud Chamber Measurements for all Experiments 
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APPENDIX C 


Raw Data Produced by the PVM for all Experiments 
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Analyzed Data 
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